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We present an analytical expression for the static many-body local field factor G-(q) of a homo- 
geneous two-dimensional electron gas, which reproduces Diffusion Monte Carlo data and embodies 
the exact asymptotic behaviors at both small and large wave number q. This allows us to also 
,_i . provide a closed- form expression for the spin-antisymmetric exchange and correlation kernel K J c (r) 

which represents a key input for spin-density functional studies of inhomogeneous electronic systems. 
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The static spin-spin response function Xs(q) of a paramagnetic electron gas (EG) can be written in terms of the 
QJ • Lindhard function xo(l) by means of the spin-antisymmetric many-body local field G-(q) through the relationship 
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where fis is the Bohr magneton. Thus G_ (q) is a fundamental quantity for the determination of many properties of a 
general electron system. It is a key input, together with the charge-charge local-field factor G+(q), in the spin-density 
functional theory (SDFT) of the inhomogeneous electron gaaJ'B and in studies.Qf quasiparticle properties (such as the 
effective mass and the effective Lande g-factor) in the electronic Fermi liquicfm 

For what concerns SDFT calculations, a common approximation to the unknown exchange-correlation energy 
functional E xc [n, £] of electron density n and spin-polarization £ appeals to its second functional derivatives, 



o 
o 



*r'(n,C;|r-r'|)- 



(2) 

n— n, C— C 



<W(r)<W(r') 

o: _ 

l/"") , where n = n-t+n i is the average local density of the EG and £ = (n^—n^/n is the average local spin-polarization. The 
C ■ local field factors G aa i (q) and the matrix of the exchange-correlation kernels are simply related in Fourier transform 
O ■ by 

"^3 ■ K! («) = / d d r e- % qr K™ (r) = -v q G aa , (q) , (3) 

a ; J 

where d is the dimensionality of the system and v q is the Fourier transform of the Coulomb potential e 2 /r. In the 
_JL , paramagnetic state the following relation holds between the charge-charge and spin-spin local-field factors G± (q) and 

£ • Gw(g): 

G ± ( g ) = GTT(g)±Gu(g) - (4) 
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Thus, from Eq. (0) and (0), one obtains that the spin-symmetric and the spin-antisymmetric exchange and correlation 
5_i ' kernels K^ c (r) are related to G± (q) by 
d | 

KtM) ^ ^M±E1M s j^e-^Kt{r) = -v q G ± (q) . (5) 



In what follows we shall only consider K~ a (r) in the case of two spatial dimensions, with d — 2 and v q = 2ire 2 /q. An 
analytical expression for K^ c (r) has already been given in Ref. |5|. _ 

A number of exact asymptotic properties of G_ (q) in two dimensions are readily proven. In particularQ, 
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where kp = \j2rrm, = y2/r s aB is the Fermi wave number, r s — yiraa| is the usual EG density parameter with clb the 
Bohr radius, \o — Ms rn / 7T ^ 2 is the (Pauli) spin susceptibility of the ideal Fermi gas while \ is the spin susceptibility 
of the interacting system. By making use of the thermodynamic definition of x we can write 
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where e c is the correlation energy per particle in the spin-polarized fluid. Unfortunately the carrelation energy at 
finite spin polarization is poorly known. By extrapolating the Diffusion Monte Carlo (DMC) datalll for G_(g) to q = 
and by observing from Eq. (@) that A_ — > 1/tt for r s — > 0, we obtain the following parametrization of A- in the range 
< r, < 10: 



A_ = 

■k + 1.4954 r s + 0.3193,/; 

The asymptotic behaviour of G_ (q) at large q is also known exactl; 

G-(q) ~ C_ -i- + S 
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for g — > oo, where C_ is proportional to the difference in kinetic energy between the interacting and the ideal gas, 
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The correlation energy per particle e c (r s ) in the paramagnetic fluid is available from Monte Carlo datacl. Finally 
B- = (7(0), g(0) being the value of the pair-correlation function at the origintEl. 

In this work we fit the values of G_ (q) originally obtained by DMC in Ref. in such a way as to obtain analytical 
expressions for both K~ (q) and K~ (r) . The functional form previously employedo for the charge-charge local- field 
factor G+(q) is also suitable for G_(g): 
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where q = q/k F , «-(r«) = ^/T~- KT0082rf and P-(q) is the polynomial P_ (q) = h 2 q 2 + h A q i + h 6 q 6 + h s q 8 . The free 
parameters a_, /?_ and /12, ..., /is ar e fitted so as to minimize the differences from the DMC numerical results. As for 
G+(q), it proves useful to have a continuous parametrization of these parameters at least in the range < r s < 10. 
We propose the following: 



a_(r 8 ) = exp 



0.2231 + 81.2115 (r s /10) 

1 + 54.6665 (r s /10) + 50.7534 (r s /10) 2 

/3_(r 8 ) = 0.8089 - 0.4025 (r s /10) 3 - 0.0941 (r s /10) 1/2 , 
, . , / 12.6262 + 20.9673 (r s /10)\ 
h2 ^- c ^{ 1 + 12.4002 (r s /10) ) ' 

h A {r s ) = 0.0531 (l - e~ 2 - 154r2 A - 0.4984 (r s /10) 3/2 + 0.4021 (r s /10) 2 , 

heirs) = -0.0076 (l - er r A + 0.0977 (r s /10) 3/2 - 0.0726 (r s /10) 2 , 
ha(r.) = -0.0027 (r./lO) . 
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In Figure [T] we compare the fit given by Eq. ( |l2|) with the DMC data for r s — 1,2,5 and 10. In Figure g we show the 
local field factor G-(q) as from Eq. (n2h for various values of r s . 



We turn next to the evaluation of K xo (r). From Eq. (|^) and (jl2|), the expression of the spin- antisymmetric exchange- 
correlation kernel in real space (in Ryd) is readily obtained as 

x(2) (r \ exp( - B-k F r/(A-K-(r s ))) , , 4 
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where JVi = -4^C_/r s , 7V 2 = -2y/2B-/r a , 7V 3 = -4V2A_ (1 - «_(r s ))/r s , N 4 = V2C-/r s f3- and AT 5 , n = 
—2a /i„ /r s . The function F n (a, x) is given by 

F n (a,x) = j o °° dvy n Jo(xy)e- a f = \ a -W r(±±^) ^ (I±^ ; l;_g) (15) 

where T(z) is Euler's Gamma function and iFi(a]b; z) is Kummer's function. In practice the function F n (a,x) can 
be obtained via the recursive relation 
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where I„(z) is the modified Bessel function of order n. It is also useful to recall that dlo(z)/dz = li(z) and that 
dli(z)/dz — Iq(z) — l\(z)/z. 

In Figure we show the matrix elements K°° (r) of the exchange-correlation kernel for various values of r s . These 
kernels describe the local structure (Pauli-Coulomb hole) of the paramagnetic electron fluid around an electron of a 
given spin. KJJ(r), being the sum of K*(r) and K~(r), has the same behavior of these two components, namely no 
structure at both small and high r s . The situation is different for Kjj (r) , which is instead given by the difference of 
K^ c (r) and K~ (r) . Although these two components are structureless at small r s , the difference between them around 
kpr ~ 1 is large and this is reflected in the large structure shown in Figure |3[ 

In conclusion, we have presented an analytic parametrization of the local field factor entering the spin response of 
the two-dimensional electron gas in the paramagnetic state, incorporating the known asymptotic behaviors and giving 
an accurate representation of tha-available quantum Monte Carlo data. We have obtained from it and from our earlier 
results on the dielectric responses analytic expressions of the exchange-correlation kernels for spin-density-functional 
calculations on inhomogeneous two-dimensional electronic systems. 
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FIG. 1. The local field factor G-(q) for various values of r s as computed according to Eq. (113), in comparison with the DMC 
data of Ref. 7 
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FIG. 2. The local field factor G-(q) as from Eq. (|12J) for various values of 
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FIG. 3. The matrix elements K]}(r) and Klj(r) of the exchange-correlation kernel for various values of r a 
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